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We discuss the analogy of the behaviour of films and drops of liquid on a rotating 

O 

\ horizontal cylinder on the one hand and substrates with regular one-dimensional 

wettability patterns on the other hand. Based on the similarity between the 
respective governing long-wave equations we show that a drop of partially wetting 
J>^| liquid on a rotating cylinder undergoes a depinning transition when the rotation 
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speed is increased. The transition occurs via a sniper bifurcation as in a recently 
described scenario for drops depinning on heterogeneous substrates. 

The manuscript was accepted by the J. Fluid Mech. in October 
2010 and will be published in the following months. 

I. INTRODUCTION 

Motivated by the question how much honey can be kept on a breakfast knife by rotating 
the knife about its long axis, Moffatt studied the destabilisation of a film on the outside of a 
rotating horizontal cylinder and the subsequent development of regular patterns of azimuthal 



nngs m . He a!so gives a lubrication equation that incorporates gravity and viseostty effects 
but neglects surface tension and inertia. It is based on a model by 40 which includes 

n 

surface tension. Related questions had been studied before by [56| employing a rotating 
lorizontal cylinder which is at the bottom immersed in a liquid bath, i.e., in contrast to 



29l | the amount of liquid on the cylinder is not an independent control parameter. The 
related case of a liquid film covering the inner wall of a horizontal rotating cylinder was 
first studied by [37| who discussed experiments performed with water employing an inviscid 
model incorporating gravitational and centrifugal forces. 

Since the early experimental results engineers and scientists alike have studied flows of 
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free-surface films of liquid on the inner or outer wall of resting or rotating horizontal cylinders 
in a number of different settings. The studied systems are of high importance for several 
coating and printing processes where rotating cylinders transport the coating material in the 
form of liquid films. In the present contribution we will discuss the analogy between film flow 
and drop motion on a rotating horizontal cylinder on the one hand and on heterogeneous 
substrates with regular wettability patterns on the other hand. First, however, we will give 
a short account of the seemingly disconnected fields. 

The theoretical analyses given by [56] and 37] for the flow on/in a rotating cylinder 
were based on linear stability considerations based on the full hydrodynamic equations for 
momentum transport, i.e., the Navier-Stokes equations. Their studies were extended later 
on, e.g., by who considered inviscid film flow on the in- and outside of the cylinder 



in a unified manner, and 
Reynolds numbers. 



45] whose analysis includes viscous film flow at small and large 



However, most ana 



yses of the non-linear behaviour are based on a long-wave or lubri- 



cation approximation 33] valid for the case where the thickness of the liquid film h(6,t) is 



small as compared to the radius of the cylinder R. The resulting evolution equation for the 
film thickness profile 



d t h = -d e {h 3 [ad e {d ee h + h) - 7cos(0)] + h} 



(1) 



was first given by 



40] and used by [29] (without capillarity effects, i.e., a = 0). In the 



dimensionless Eq. ([TJ 



a 



3e 3 a/Rojfj and 7 = e 2 pqR/3wf] are the scaled dimensionless 



18] . It is important to note that both 



surface-tension and gravity parameter, respectively 
depend on the angular velocity of the rotation u, the dynamic viscosity fj, the cylinder 
radius R, and the ratio e of mean film thickness h and R. The remaining material constants 
are surface tension a and density p, and g is the gravitational acceleration. The angle 9 
determines the position on the cylinder surface and is measured anti-clockwise starting at 
the horizontal position on the right. The force 7cos(0) corresponds to the component of 
gravity parallel to the cylinder surface. 



The 



e.g. 



18 



ong-wave equation (CQ) including surface tension effects has been studied in detail by, 



41] and 23] . Note that for the parameter ranges where Eq. (0Q) applies, it is va 



without difference for films on the outside and inside of the cylinder (cf., e 



. g ., 0,Q, 



id 



ail). 



A preliminary comparison of the long-wave approach to results obtained with the Stokes 
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equation is given by [36]. Various extensions of Eq. §T§ were developed. See, in particular, 
[l| for a discussion of hydrostatic effects. Other higher order terms related to gravity, inertial 



and centrifugal effects have been included 
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25. 
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32] . In particular, [3l| and 



25] present a systematic derivation of several such higher order models and give also detailed 
comparisons to a number of previous studies as e.g., {32]] , Q and fl. The film stability 
and non-linear evolution of a film on the outside of a non-isothermal horizontal cylinder was 



considered without 



43] and with 13] rotation. Most of the aforementioned studies consider 



two-dimensional situations. The fully three-dimensionional problem is studied for the cases 



without 



551 ] and with 
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iquid drops and films on heteroge- 



30] rotation. 

The second system that we would like to discuss are 
neous solid substrates as studied experimentally, e.g. by [46] and [42]. Although in principle 
the heterogeneity may result from substrate topography or chemical heterogeneities we focus 
here on the case of a smooth flat substrate with chemical heterogeneities, i.e., we assume 
the substrate wettability depends on position. Horizontal heterogeneous substrates are often 
considered in connection with micro-patterning of soft matter films via thin film dewetting 



Sin 
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44|. Another application are free-surface liquid channels in microfluidics 



16]. 



Often one employs thin film evolution equations obtained through a long-wave approxima- 
tion to study the dynamics of dewetting for an initially flat film or to analyse steady drop 
and ridge solutions and their stability on such patterned substrates (see, e.g., and 49]). 
Static structures may also be studied using a variational approach 10, 127 ]. 

More recently a model system was introduced to study the dynamics of drops on het- 
erogeneous inclined substrates. The models apply also to such drops under other driving 
forces along the substrate 52]. The corresponding dimensionless evolution equation for the 
one-dimensional film thickness profile (i.e., describing a two-dimensional drop) is of the form 



d t h= -d x {h 3 [d x (d xx h + U(h,x))}+jl} 



(2) 



where fl(h) represents the lateral driving force (fih 3 in the case of an inclined substrate, where 
/i is the inclination angle). The position-dependent disjoining pressure U(h,x) models the 
heterogeneous wettability of the substrate. One notices at once that Eqs. (|2]) and (CQ) are 
rather similar: Identifying x = 8, £i(h) = h and U(h,x) = ah — 7sin(x), Eq. ([[]) may for 
a = 1 be seen as a special case of Eq. ((2]). Here, we use the analogy to study drops or films 
on a rotating cylinder. In particular, we show that drops on a rotating cylinder may show 
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rather involved depinning dynamics analogeous to the behaviour of drops on heterogeneous 
substrates ol. I51L 52|. 

To illustrate the behaviour we choose the case of drops of a partially wetting liquid on 
a substrate with a periodic array of localised hydrophobic wettability defects [52]. On a 
horizontal substrate one finds steady drops that are positioned exactly between adjacent 
defects, i.e., the drops stay on the most wettable part. At a small lateral driving force \x the 
x —x symmetry of the system is broken. However, the drops do not slide continuously 
along the substrate as they do in the case of a substrate without defects 53] . Their advancing 
contact line is blocked by the next hydrophobic defect, i.e., one finds steady pinned drops 
at an upstream position close to the defect. As the driving force is increased, the drops are 
pressed further against the defect. As a result the steady drops steepen. At a finite critical 
value n c the driving force allows the drop to overcome the adverse wettability gradient 
and the drop depins from the defect and begins to slide. In a large part of the parameter 
space spanned by the wettability contrast and the drop size the point of depinning at /i c 
corresponds to a Saddle-Node Infinite PERiod (sniper) bifurcation. Beyond fi c the drops 
do not slide continuously down the incline. They rather perform a stick-slip motion from 



defect to defect 



51 



521 ] . However, in another region of the parameter space the depinning 



can occur via a Hopf bifurcation, i.e., with a finite frequency. In particular, this happens in 
the case of rather thick wetting layers. 

The aim of the present contribution is to elucidate under which conditions a similarly 
intricate 'depinning' behaviour can be found for a pendant drop underneath a horizontal 
cylinder when increasing the speed of rotation. In the following section [III we will discuss 
our scaling and introduce a long-wave model for a film of partially wetting liquid on a 
rotating cylinder. It allows (even without gravitation) for the coexistence of drops and a 
wetting layer. We discuss several physical situations it can be applied to. Next, section 
IIHI presents steady state solutions depending on the non-dimensional rotation velocity for 
several thicknesses of the wetting layer and discusses the depinning scenario in the case of 
the partially wetting liquid. Finally, we conclude and give an outlook in Section [TV] 
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II. MODEL 

Before we introduce our long-wave model we discuss our scales. In the non-dimensionali- 
sation used in most works studying Eq. (PQ) the angular velocity of the cylinder rotation u 
is not reflected in a single dimensionless parameter but both dimensionless parameters are 



inversely proportional to u. Furthermore, the time scale depends on u [18|, |40|. This and 
a further re-scaling of the equation for steady drop and film states, that absorbs the flow 
rate into both dimensionless parameters and also into the thickness scale, implies that the 
existence of two steady solutions for identical parameters (thin capillary film and pendant 



drop with through-flow), shown analytically by 4]j and numerically by [23], actually refers 
to solutions of different liquid volume. 

Such a scaling is not convenient for our study as we would like to investigate the change in 
system behaviour for fixed liquid volume when increasing the rotation speed of the cylinder 
from zero, i.e., u should only enter a single dimensionless parameter that will be our main 
control parameter. A scaling fit for our purpose is the one based on the scales of the gravity- 
driven drainage flow employed by 



Note also that in particular [23[ and to a lesser extent 55] represent some of their solutions 
in a way that misrepresents their findings and, in general, the concept of an equation in long- 
wave approximation. Whenever one shows the solutions as a thickness profile on a cylinder 
one has to pick a particular radius R and smallness ratio e for the representation (cf. Figs. |3j 
|5] and [9] below). This is a rather arbitrary choice as the equation in long- wave approximation 
is strictly valid only in the limit e — > 0, ie. it is always only an approximation to a real 
physical situation where e is finite. 

In the present work we focus on solutions without axial variation, i.e., we consider the 



two-dimensional physical situation depicted in Fig. [TJ Following lj| we use h, R, pgh 2 /rj, 
epgh 2 /i], 3r]/ehpg and pgh to scale the coordinate z orthogonal to the cylinder surface, 
the coordinate x = 9R along the surface, the two velocity components v x and v z , time, and 
pressure, respectively. The ratio of z— and x— scale is the smallness parameter e = h/R, i.e., 
the film height is scaled by h — eR. Here the angle 9 is defined in clockwise direction starting 
at the upper vertical position, i.e., a pendant droplet underneath the resting horizontal 
cylinder has its centre of mass at 9 = ir. 

The evolution equation is derived from the Navier-Stokes equations with no-slip boundary 
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FIG. 1: Sketch of a drop of partially wetting liquid coexisting with a thick wetting layer on a 
rotating cylinder. 

conditions at the (rotating) cylinder, and tangent- and normal force equilibria at the free 



surface using the long-wave approximation 



14 



33 



40 



431 ] . Note that here we furthermore 



assume that the ratio h/R and the physical equilibrium contact angle are both of 0(e) 
(cf. Section HVl) . This gives 



ct h 



d ft lh 3 de 



Bo" 1 (d ee h + h) - cos(fl) + 11(h) + fi/ij , 



where r is the non-dimensional time, and 

Bo fl3 « 



and Q 



r/uR 



> 



(3) 



(4) 



ha ph 2 g 

are an effective Bond number and a rotation number for a clock-wise rotation, respectively. 
The latter corresponds to the ratio of the rotation velocity at the cylinder surface and the 
drainage velocity. The function 11(h) is a non-dimensional disjoining pressure that accounts 



for the wettability of the liquid film on the cylinder surface js, |l7|, |21 |. 

To make the analogy with the depinning droplet on a heterogeneous substrate more 
obvious we introduce another timescale t = r/Bo, disjoining pressure II = BoLT and rotation 
number Q = BoQ and obtain 



t h = -d e {h 3 d e [d ee h + h- Bocos(fl) + 11(h)] + tth} 



with 



n 



rjuR 
e 3 a ' 



(5) 



(6) 
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Defining the derivative 

d h f = -Tl(h) -h + Bo cos(6) (7) 

of a 'local free energy' f(h,6) and identifying fx = Qh, Eq. (jSJ) takes the form of Eq. (J2]) 
and the Bond number takes the role of a 'heterogeneity strength'. Note that the scaling 
used is well adapted for discussing independently the influences of rotation and gravity in 
the presence of capillarity. However, it does not allow for a study of the limit of vanishing 
surface tension. 

As disjoining pressure we choose here the combination of a long- and a short-range power 
law U(h) = Ha//i 3 (1 — b/h 3 ), where Ha is a dimensionless Hamacker constant. For b > 
one has a precursor film or wetting layer thickness ho = b 1 ^ 3 38|, |39|, |48| . This is valid only 
for Ha< 0, i.e., when the long-range interaction is destabilising. This corresponds to the 
case of a partially wetting liquid. For Ha> 0, a rather thick wetting layer can be stabilised 
and the b/h 6 term is of no further relevance and could just as well be dropped. 

Note that dhf may contain other terms beside the ones in Eq. (JTj). For a heated or cooled 
cylinder a term (3/2) BoBi Ma log[/i/(l + Bi h)\ + 1/(1 + Bi h)\ would enter where Bi is the 



Biot number and Ma an effective Marangoni number |34j . 1431 . [50[. If the rotating cylinder 
forms the inner electrode of a cylindrical capacitor, a term Vo/[h + (d — h)e} 2 has to be 
included for a dielectric liquid and DC voltage of non-dimensional strength Vo (e is the 
electric permittivity and d the distance between the two cylinders (see appendix of 52| and 



22]). The similarity in the behaviour of depinning drops under the influence of different 



ohysical effects in the case that the effective pressure terms 'look similar' is discussed by 



52 1 . 



Although here we do not consider such other physical effects explicitly, we will explore 
the effect of corresponding thick wetting layers by adapting the parameters of our disjoining 
pressure accordingly. A thick wetting layer that coexists with droplets might, for instance, 
result from the interplay of a destabilising thermal or electrical effect and a stabilising long- 
range van der Waals interaction. We obtain adequate values of Ha and b for our present 
n(/i) by relating them to the (non-dimensional) thickness of the wetting layer ho and a static 
macroscopic contact angle f3 by [see, e.g.. Il7| 

b = h 3 and Ha = -^h 2 . (8) 

Note, that (3q is the angle in long- wave scaling, i.e., a small physical equilibrium contact 
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angle (3 eq = e(3 corresponds to a long- wave contact angle (3 of 0(1). In the limit of a non- 
rotating cylinder (i.e. Q = 0) and assuming the disjoining pressure corresponds to complete 
wetting (Ha> 0, 1/h 6 term may be dropped), our equation ([5]) corresponds to the one given 



by |43j for the isothermal case. There, however, first a viscous scaling is used, which they 
transform in a second step - a rescaling of time and therefore velocities with the Galileo 
number - into the 'drainage scaling' used here. Our equation corresponds also to that of 



181 ] when the disjoining pressure is added to their equation. 

To analyse the system behaviour we determine in the following steady-state solutions 
that correspond, e.g., to pendant droplets. The behaviour beyond the depinning thresh- 
old is analysed using time-stepping algorithms. As we restrict ourselves to the two- 
dimensional physical situation an explicit scheme for stiff equations suffices for the latter. 
The steady-state solutions are obtained using the continuation techniques of the package 



AUTO 



121 ] . The steady and time-periodic solutions are characterized by their L 2 norm 



\\Sh\ \ = \J (l/27r) Jq 71 \h{6) — l) 2 dO and time- averaged L 2 norm 



\\5h\ \ = yj (1/2ttT) Jq J q 27T \h{9) — l) 2 d8dt, respectively. T is the time period. 

Note finally that Eq. ([5]) is the result of our choices for the relative order of magnitude of 
the involved dimensionless numbers that is based on the choice of physical effects that shall 



be discussed (cf. |33|). Once this is accepted, equation fl5]) is to 0(1) asymptotically correct. 
Its form as a conservation law implies that J Q 27r hdO = 2n. To 0(1) this corres pon ds to mass 
conservation. For a discussion of higher order corrections to this picture see |25| . 



III. RESULTS 



A. Parameters 



We base our estimate for realistic Bond and rotation numbers on two liquids typically 



studied in the literature: (i) water at 25°C as in 43] and a silicone oil with o = 0.021N/m, 



1] = lkg/ms and p = 1200kg/m 3 as in 14j. As experimentally feasible configurations we 
assume angular rotation velocities u> = 0.1 .. . 10s" 1 , cylinder of radii R = 10~ 3 . . . 10~ 2 m 
and smallness ratios e = 0.01... 0.1. It is only for the purpose of illustration of steady 
droplets on a cylinder that we later use e — 0.1 in selected figures. 

For oj = Is" 1 and a smallness ratio e = 0.1 we find for a cylinder of radius R = 10 _2 m for 
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material (i) Bo= 136 and Q = 0.14, whereas material (ii) gives Bo= 561 and Q = 476. For a 
smaller cylinder radius of R = 10" 3 m one has (i) Bo= 1.36 and Q = 0.014, and (ii) Bo= 5.61 
and Q = 47.6. Decreasing the smallness ratio to e = 10 -2 increases all Bond numbers by 
a factor 10 and all rotation numbers by a factor 10 3 . Focusing on cylinder diameters in 
the millimetre range, we mainly investigate Bo< 10. Note that experiments may also be 
performed with larger cylinders. To keep the Bond number in the interesting range one could 
decrease the relevant density by replacing the ambient gas by a second (immiscible) liquid. 
A feasible experimental set-up could be a horizontal Taylor-Couette apparatus with rotating 
inner cylinder. The theoretical framework would also need to be amended - replacing the 
one-layer theory [Eq. (jHJ)] by a closed two-layer model similar to 28^ . 

Note, that silicone oils may have viscosities 2 orders of magnitude smaller or larger than 
the chosen value and one is able to vary the velocity of rotation in a wide range. This 
implies that there is some flexibility in the choice of Q. We will find that for Bo= 0(1) the 
interesting range for the rotation number is Q = 0(1). For instance, with Bo= 1 one finds 
that depinning occurs at Q = 1.68. 

As the present work aims at establishing the qualitative correspondence between the 
behaviour of a depinning drop on a heterogeneous substrate and the one of a droplet on a 
rotating cylinder we restrict our results to one partially wetting case (choosing ho — 0.1 and 
Po = 2) and the completely wetting case (/3 = 0). For comparison, selected results are also 
shown for the contact angle /3q = 1. 

Note that here the scaling and therefore the discussion of the smallness parameter e is 
based on e = h/R, i.e., we use the ratio of mean film thickness and cylinder radius. In the 
case of pendant drop solutions one also has to keep the maximal film height h max i n mind 
when discussing e as the scaling becomes questionable when /i max ^> h. This is, however, 
normally not the case. In the present work the ratio h m!ix /h is always below 3. Although 
an inspection of the figures in 23[ shows a ratio h max /h that is not much larger, one notices 
that there e itself (based on h) seems to be larger than one (see remark above Eq. (JIB ). 
The same applies to Fig. 3 of 55]. There exists no such problem in the region where only 
a thin wetting layer of thickness h min covers the cylinder even if h min <C h. Actually, the 
long-wave approximation gets better there. In an approach based on matched asymptotics 
one could further simplify the governing equation in this region (cf. (7|). However, such an 
approach is not taken here as it is of limited use when considering parameter regions that 
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FIG. 2: Solutions on a horizontal cylinder without rotation (Q = 0) are characterised by their 
L 2 norm as a function of the Bond number for two different equilibrium contact angles (3q and 
ho = 0.1, h = 1.0. 

show qualitative changes in solution behaviour, e.g., close to the sniper bifurcation discussed 
below. 



First, we consider the case of a horizontal cylinder without rotation (Q = 0) and determine 
steady drop and film solutions. Employing the Bond number as control parameter we 
obtain families of steady state solutions for two selected values of the equilibrium contact 
angle /3q. Inspecting Fig. |2]one notices that above a critical value of the Bond number Bo c 
there exists only a single solution. It corresponds to a symmetric pendant drop located 
underneath the cylinder. An example for Bo= 1 and /3q = 2.0 is given as solid line in Fig. |3j 
However, two additional steady solutions exist below Bo c . One of them corresponds to an 
unstable symmetric drop sitting on top of the cylinder (dotted line in Fig. [3]). The other 
one corresponds to an unstable solution that has two minima - a deep one underneath the 
cylinder and a shallow one on top of it (dashed line in Fig. |3]). 

We find that Bo c becomes smaller with decreasing contact angle /3q, he., the range of Bo 



B. Partially wetting case 
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FIG. 3: Steady film thickness profiles on a horizontal cylinder without rotation (f2 = 0) for Bo= 1 
and fio = 2.0 (corresponding to marks (i) to (iii) in Fig. [2]). The remaining parameters are as in 
Fig. [2j Panel (a) gives the profiles as h(6), (b) gives the space-dependent part of dhf (Eq. 0), 
and panel (c) shows the profiles on a cylinder. Note that for illustration purposes we have (rather 
arbitrarily) assumed a radius R = 10, i.e. e = 0.1 (cf. discussion in Section [III Aj) . Line styles in 
(c) correspond to those in (a). The cylinder surface is represented by the solid black line. 

numbers where steady structures beside the pendant drop exist becomes smaller. The com- 
pletely wetting case (/3 = 0) is qualitatively different and is discussed below in Section UlI CI 
Increasing (3 above (3 = 2 the critical value Bo c strongly increases, e.g., for f3 = 10 one 
finds Bo c = 43.3. 

Note, that the pendant drop underneath the cylinder [drop on top of cylinder] corresponds 
on a horizontal heterogeneous substrate to the drop on the more [less] wettable region as 
discussed by . The limit of zero Bond number is analogue to the horizontal homogeneous 
substrate. 

Next, we increase the dimensionless angular velocity Q for several selected Bond numbers 
and determine the steady state solutions (Fig. H]). Focusing first on Bo= 1 we recognise at 
Q = three solutions as marked by the vertical dotted line in Fig. El The curve of largest 
norm corresponds to the pendant drop. As Q increases from zero, its norm and amplitude 
decrease slightly, whereas its centre of mass moves towards larger 9, it is 'dragged along' by 
the rotation without changing its shape much. This can be well appreciated in Fig. [5] which 
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FIG. 4: Solutions on a horizontal rotating cylinder as a function of the Rotation number £1 for 
various Bond number Bo as given in the legend and contact angle /3o = 2.0. The drop and film 
profiles are characterised by their I? norm. The remaining parameters are ho = 0.1, and h = 1.0. 




FIG. 5: Steady film thickness profiles on a rotating horizontal cylinder for Bond number Bo= 1 
and /3q = 2.0 for Rotation numbers O as given in the legend. All shown solutions are situated on 
the stable upper branch in Fig. HI The remaining parameters are as in Fig. HJ Panel (a) gives the 
profiles as h(9), (b) gives the space-dependent part of dhf (Eq. (J7J)), and panel (c) illustrates the 
profiles on a cylinder assuming a radius R = 10. Line styles in (c) correspond to those in (a). 
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FIG. 6: Bifurcation diagram for the depinning transition for a drop on a rotating cylinder with 
Bo= 1 and f3o = 2.0. Shown are the L 2 norm for the branch of steady solutions as obtained by 
continuation (solid line), selected steady solutions as obtained by direct integration in time (circles) 
and the time-averaged 1? norm for the unsteady solutions beyond the sniper depinning bifurcation 
(triangles). The inset gives for the latter a log- log plot of the dependence of the inverse temporal 
period on the distance from the bifurcation fi — J7 sn . The solid line indicates an exponent of 1/2. 

gives selected profiles for the branch of pendant drop solutions. The branch terminates in a 
saddle-node bifurcation at Q sn = 1.68 where it annihilates with one of the unstable branches. 
The third branch, i.e., the one of lowest norm, shows for increasing Q an increase in the 
norm, then it undergoes two saddle-node bifurcations at — 1.313 and Q = 1.329. At 
larger Q the norm decreases again. Note that the two saddle-node bifurcations are a first 
sign of rather complicated re- connect ions which occur when one increases Bo from 1.0 to 
2.5. The re-connections involve further (unstable) branches that are not shown in Fig. [Hand 
result in the loop structure observed for Bo= 2.5 at about Q = 1.5. A further indication for 
the existence of additional branches is the observation that the value of the norm at Q = 
for Bo= 1.0 in Fig. H] does not agree with the lowest norm at Bo= 1.0 in Fig. |2j As the 
re-connections are not relevant for the depinning transition we will here not consider them 
further. 

At Bond numbers smaller than Bo= 1 the behaviour is qualitatively the same as described 
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(a) 



h 2. 




(b) 



2.69 



FIG. 7: Space-time plots illustrating the time evolution of co-rotating drops beyond depinning via 
a sniper bifurcation (at O > sn ). Shown is one period in space and time (a) close to depinning 
at Rotation number 0, = 1.8 with a temporal period of T = 14.8, and (b) far from depinning at 
= 3.0 with T = 2.7. The Bond number is Bo= 1 and j3 = 2.0. 

for Bo= 1. However, the critical Q sn becomes smaller with decreasing Bo. Above, we have 
discussed that the behaviour at Q = changes qualitatively when decreasing Bo beyond 
Bo c (Fig. [2]). Although this also affects the behaviour at small Q it has less influence on 
the behaviour at larger fl (Fig. HJ). In particular, the saddle-node at Q sn where the pendant 
drop solution ceases to exist, persists for larger Bond number as does the related dynamic 
behaviour (see below). Note that for Bo= 10 the saddle-node is at Q sn « 6.2 outside the 
range of Fig. HJ 

An interesting resulting question is what happens to a stable obliquely pendant droplet 
at Q < Q sn when the rotation velocity is further increased such that it is slightly larger than 
n sn ? Does the solution approach the remaining steady solution? Simulations of Eq. (J5]) 
show that this is not the case. For Q > Q sn the droplet moves in a non- stationary way with 
the rotating cylinder and corresponds to a space- and time-periodic solution. The resulting 
bifurcation diagram for Bo= 1 is given in Fig. [6j It shows that the 'new' branch of space- and 
time-periodic solutions (characterised by its time-averaged norm) emerges from the saddle 
node bifurcation of the steady solutions at Q sn . The behaviour close to and far from the 
saddle-node bifurcation is illustrated in the space-time plots of Fig. [7J Panel (b) shows 
that far away from the bifurcation the drop moves continuously with its velocity and shape 
varying smoothly. It moves fastest when its maximum passes 9 w 7r/2 (velocity w 4.6) 
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and slows down when passing 9 rs 3ir/2 (velocity ps 1.2), i.e., respectively, when gravity 
most strongly supports and hinders the motion driven by the cylinder rotation. The ratio 
of largest to smallest velocity is about 4:1. When approaching the bifurcation from above, 
the droplet still moves with the rotation but the time scales of the slow and the fast phase 
become very different. For instance, at Q = 1.8 (Fig. [7] (a)) the drop moves fastest when 
passing 9 « 7r/2 (velocity ~ 3.3). However, when the drop is situated at about 9 = 37r/2 
it barely moves (velocity w 0.1) resulting in a velocity ratio of 33:1 between the fastest 
and slowest phase. The resulting overall behaviour strongly resembles the stick-slip motion 
discussed in the context of contact line motion on heterogeneous substrates j^. Note that 
it also resembles so-called sloshing modes found for a partially liquid-filled rotating cylinder 
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54j. 



The ratio of the velocities diverges when approaching the bifurcation point. In conse- 
quence, the frequency XjT of the periodic drop motion goes to zero. As shown in the inset of 
Fig. El the dependence corresponds to a power law 1/T ~ (O — fisn) 1 / 2 . This indicates that 
the bifurcation at fl SQ is actually a Saddle-Node- Infinite-PERiod bifurcation (or 'sniper' for 
short; see J^, and discussion by 52]). 

The described behaviour is quite generic for the partially wetting case, i.e., it is found 
in a wide range of parameters, in particular, Bond number, contact angle and wetting 
layer thickness. However, the behaviour changes dramatically when strongly decreasing the 
contact angle, i.e., when approaching the completely wetting case. 



C. The completely wetting case 

After having discussed the intricate depinning behaviour in the partially wetting case, 
to emphasise the contrast, we consider briefly the completely wetting case (/3o = 0). Note 
that the physical setting is then identical to the one employed in most studies of the classic 



Moffatt problem 
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40j. However, the parametrisation used here is different and results 
in Eq. ([3]) with Ha= 0. We analyse the system behaviour as above by studying the steady 
state solutions without and with rotation. 

In the case without rotation, the first difference one notices is that for fixed Bond number 
only one solution exists. It corresponds to the pendant drop solution. With increasing Bo, 
i.e., with increasing importance of gravity, its amplitude monotonically increases. Fixing 
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FIG. 8: Solutions for a film of completely wetting fluid (/3q = 0) on a horizontal rotating cylinder 
as a function of the Rotation number for various Bond numbers Bo as given in the legend. The 
solutions are characterised by their L 2 norm. The remaining parameter is h = 1.0. 




FIG. 9: Steady film thickness profiles for the case of complete wetting with Bo= 1 for Rotation 
numbers f2 as given in the legend. The remaining parameters are as in Fig. [8l Panel (a) gives the 
profiles as h(9), (b) gives the space-dependent part of dhf (Eq. ((7J)), and panel (c) illustrates the 
profiles on a cylinder assuming a radius R = 10. Line styles in (c) correspond to those in (a). 
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Bo and increasing the rotation velocity from zero, one finds a monotonic decrease of the 
norm (Fig. [8]) and height of the drop. Its maximum moves towards larger angular position 
6, i.e., the liquid is dragged along with the rotation. However, in the wetting case there 
is no force besides gravity that can favour drops as compared to a flat film. As gravity 
acts downwards it only favours pendant drops. Therefore, the drop does not survive as a 
'coherent structure' when dragged upwards with the rotation: it is flattened and smeared 
out. Typical profiles are given in Fig. |9j Note that we do not show the exact case without 
rotation, but choose Q = 1CT 10 as then a dynamically created wetting layer still exists that 
is numerically advantageous (cf. [50]). 

When increasing Q further than shown in Fig. E] one finds that the thickness profiles 
approaches a flat film. For instance, for Bo= 1 the norm decreases below 10~ 2 at Q ~ 70. 

IV. CONCLUSION 

Based on the observation that the equations in long-wave approximation that govern film 
flow and drop motion (i) on or in a rotating cylinder and (ii) on a heterogeneous substrate are 
rather similar, we have explored whether the analogy can be exploited, i.e., whether it allows 
results obtained for one system to be transferred to the other one. In particular, we have 
found that indeed for drops on a rotating cylinder there exists a counterpart of the rather 
involved depinning dynamics described recently for drops on heterogeneous substrates. 

To study the effect we have introduced an alternative scaling. This has been necessary 
because the commonly used scaling contains the angular velocity of the cylinder rotation in 
both dimensionless parameters and in the time scale. Together with a further re-scaling of the 



steady state equation involving the flow rate 
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40( this does not allow for an investigation 



of either the system behaviour when increasing the rotation speed of the cylinder from 
zero or the existence of multiple solutions for fixed rotation speed and liquid volume. The 
scaling employed here does allow for such studies as the rotation speed only enters a single 
dimensionless parameter. Furthermore, it allows us to discuss the analogy between the two 
systems of interest in a rather natural way. In particular, downward gravitation and rotation 
speed for the drop on the rotating cylinder correspond to the heterogeneous wettability and 
lateral driving force, respectively, for drops on heterogeneous substrates. 

Guided by this analogy we have studied the behaviour of drops of partially wetting liquid 
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on the rotating cylinder. As a result it has been shown, how 'switching on' gravity (in- 
creasing the Bond number) without rotation changes the solution behaviour dramatically as 
the 'gravitational heterogeneity' along the cylinder surface effectively suppresses multidrop 
solutions until only the single pendant droplet underneath the cylinder survives. We have 
furthermore found that increasing the rotation from zero to small values, the drop is dragged 
along to a stable equilibrium position where downwards gravity and upwards drag compen- 
sate. At a critical speed, however, the retaining downwards force is too small and the drop 
is dragged above the left horizontal position. It then continuously moves with the rotating 
cylinder in a non-constant manner. Close to the transition the drop shows stick-slip motion. 
An analysis of the time-dependent behaviour has shown that the frequency related to the 
periodic motion goes to zero at the threshold following a power law with power 1/2. This 
and the related bifurcation diagram show that the observed transition corresponds to depin- 
ning via a sniper bifurcation. The behaviour found is rather generic for the case of partial 
wetting. Although, here we have presented results for particular parameter values, it can 
be observed in a wide range of parameters, in particular, Bond number, equilibrium contact 
angle and wetting layer thickness. However, we have also found that the behaviour changes 
dramatically when strongly decreasing the equilibrium contact angle, i.e., when approaching 
the completely wetting case that is normally studied in the literature. When increasing the 
velocity of rotation in this case we have seen a transition from pendant drops underneath 
the cylinder to a nearly uniform film around the cylinder as fl^j . 

The presented results indicate that it may be very fruitful to further explore the analogy 
of films/drops on rotating cylinders and heterogeneous substrates. We expect that the 
recent exploration of the three-dimensional case for depinning drops [5| and for the various 
transitions between ridge, drop and rivulet states |4| based on tools for the continuation of 
steady and stationary solutions of the two-dimensional thin film equation (physically three- 
dimensional system) [6[ is also relevant for the case of a rotating cylinder. A future study 
could, e.g., elucidate the relation between the formation of azimuthal rings, pendant ridges 
and sets of drops. 

The present circular cylinder is an analogy to a sinusoidal wettability profile. Cylinders 
with other cross sections, for instance, an elliptical cylinder 20|] would correspond to more lo- 
calised defects of the corresponding heterogeneous substrate. However, the analogue system 
turns out to be more complicated as a rotating non-circular cylinder introduces the aspect 
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of a time-periodic non-harmonic forcing parallel to the substrate into the heterogeneous 
substrate system. 

Our equation (jHJ) is derived under two conditions: (i) the mean (and also the maximal) 
film thickness has to be small as compared to the radius of the cylinder; and (ii) the local 
surface slope (e.g., the physical equilibrium contact angle) has to be small. Here we have 
assumed that the two related smallness parameters are of the same order of magnitude. In 
a next step one may introduce two different smallness parameters and discuss a number of 
distinct limits in dependence of their ratio. Such a systematic asymptotic study would also 
permit to discuss the influence of other effects on the depinnin g be haviour like, for instance, 



the effects of the hydrostatic pressure 



03, 



centrifugal forces 14j and inertia 19l. |25|. 
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